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Rotational bands have been observed to emerge in ab initio no-core configu¬ 
ration interaction (NCCI) calculations for p-shell nuclei, as evidenced by rotational 
patterns for excitation energies, electromagnetic moments, and electromagnetic transi¬ 
tions. We investigate the ab initio emergence of nuclear rotation in the Be isotopes, 
focusing on 9 Be for illustration, and make use of basis extrapolation methods to obtain 
ab initio predictions of rotational band parameters for comparison with experiment. 
We find robust signatures for rotational motion, which reproduce both qualitative and 
quantitative features of the experimentally observed bands. 
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1. INTRODUCTION 

Nuclei exhibit a wealth of collective phenomena, including clustering, rotation, 
and pairing [1], which we may now seek to understand through ab initio approaches. 
The challenge of ab initio nuclear theory is to quantitatively predict the complex and 
highly-correlated behavior of the nuclear many-body system, starting from the under¬ 
lying internucleon interactions. Significant progress has been made in the ab initio 
description of light nuclei through large-scale calculations, and signatures of collec¬ 
tive phenomena have now been obtained in ab initio calculations [2-6]. The emer¬ 
gence of rotational bands, in particular, has recently been observed [7-10] in ab initio 
no-core configuration interaction (NCCI) [11] calculations for p-shell nuclei. Rota¬ 
tional patterns are found in the calculated level energies, electromagnetic moments, 
and electromagnetic transitions. 

Yet, NCCI calculations are, of necessity, carried out in a finite, truncated space. 
Computational restrictions limit the extent to which converged calculations can be 
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obtained to identify collective phenomena. 

Therefore, natural questions surrounding the emergence of rotation in ab initio 
calculations include: How recognizable is the rotation (from the calculated observ¬ 
ables)? How robust is the rotation (in incompletely converged calculations)? How 
realistic is the rotation (when quantitatively compared to experiment)? How does the 
rotation arise (or what is the intrinsic physical structure)? In this proceedings contri¬ 
bution, we briefly introduce ab initio NCCI calculations (Sec. 2) and the signatures 
of nuclear rotation (Sec. 3). We then consider the emergence of rotation in NCCI 
calculations (Sec. 4). We introduce and augment the central ideas and discussions 
of Ref. [10], touching on the above questions, using 9 Be for illustration. We then 
compare the rotational energy parameters of the calculated bands in 7-12 Be against 
those of their experimental counterparts. We apply an exponential basis extrapola¬ 
tion scheme for the energies, to estimate the band parameters which would be found 
in an untruncated calculation. 


2. NCCI CALCULATIONS 

In NCCI calculations, the nuclear many-body Schrodinger equation is formu¬ 
lated as a Hamiltonian matrix eigenproblem. The Hamiltonian is represented with 
respect to a basis of antisymmetrized products of single-particle states, convention¬ 
ally, harmonic oscillator states. The problem is then solved for the full system of 
A nucleons, /.<?., with no inert core. However, calculations must be carried out in a 
finite-dimensional space, commonly obtained by truncating the basis to a maximum 
allowed number 7V max of oscillator excitations. Convergence toward the exact re¬ 
sults — as would be achieved in the full, infinite-dimensional space — is obtained 
with increasing N max . However, the basis size grows combinatorially with JV max , so 
the maximum accessible iV max is severely limited by computational restrictions. 

The calculated eigenvalues and wave functions, and thus observables, depend 
both upon the basis truncation iV max and on the length parameter for the basis, which 
is specified by the oscillator energy Ml. These dependences are illustrated for the 
calculated ground-state energy of 9 Be in Fig. 1(a), for 2 < N max < 10. For a fixed 
lV max , a minimum in the calculated energy is obtained in the range Ml ~ 20MeV to 
25 MeV, providing a variational upper bound. Then, as JV max is increased, a lower 
calculated ground state energy is obtained at each Ml. The approach to convergence 
is marked by approximate A r max independence (a compression of successive energy 
curves) and Ml independence (a flattening of each curve around its minimum). While 
a high level of convergence (at the keV scale) may be obtained in the lightest nuclei 
(in particular, see Fig. 3 of Ref. [12] for 4 He), the situation is challenging for the Be 
isotopes. Although the decrease in the variational minimum energy for 9 Be becomes 
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Fig. 1 - The fV max and Kl dependence of calculated observables for 9 Be: (a) the energy of the 3/2“ 
ground state and (b) the yrast 5/2“ —> 3/2“ electric quadrupole transition strength. Calculated values 
are shown as functions of hQ. for N max = 2 to 10 (as labeled). 


smaller with each step in iV max [Fig. 1(a)], even at iV max = 10 these changes are still 
at the MeV scale. 

For electric quadrupole transition strengths, convergence is even more elusive, 
as shown for the transition between the lowest 5/2“ state and 3/2~ ground state 
of 9 Be in Fig. 1(b). Convergence — as manifested in lV max independence and M> 
independence — would here be reflected in a compression of successive lV max curves 
and a “shoulder” in the plot of the B(E2) against hVL (see Fig. 3 of Ref. [12] for 
analogous convergence, albeit of a different observable, in 4 He). Fiowever, at most 
hints of such features are apparent at N max = 10 in Fig. 1(b). Consequently, there is 
no obvious way to extract a value for the E 2 strength (other than to arbitrarily choose 
an hQ value at which to do so, as we shall consider in Sec. 4). 

3. COLLECTIVE NUCLEAR ROTATION 

Nuclear rotation [1] arises when there is an adiabatic separation of a rotational 
degree of freedom from the remaining internal degrees of freedom of the nucleus. A 
rotational state factorizes into an intrinsic state \4 >k) and a rotational wave function 
of the Euler angles $, describing the collective rotational motion of this intrinsic state. 
(Specifically, we consider an axially symmetric intrinsic state, with definite angular 
momentum projection K along the intrinsic symmetry axis.) The full nuclear - state 
ICj/cm}■ with total angular - momentum J and projection M, has the form 

HjKm) OC jdd [ + (- ) J+K @ J M- K (mK,0)] , (1) 

Rotational Intrinsic 

where \(J)k] ’0) represents the intrinsic state \4>k) after rotation by i), and the second 
term involving the ^-conjugate state | cpf >; i)) arises from discrete rotational symme- 
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try considerations. 

The principal signatures of rotational structure reside in rotational patterns in 
the energies and electromagnetic multipole observables for the set of rotational states. 
Rotational band members, sharing the same intrinsic state but differing in J and 
hence their rotational wave functions, have energies following the rotational formula 
E(J ) = E{) + AJ(J + 1), where the rotational energy constant A is inversely related 
to the moment of inertia of the intrinsic state. However, for K = 1/2 bands, the 
Coriolis contribution to the kinetic energy significantly modifies this pattern, yielding 
an energy staggering 

E(J) = E 0 + A[J{J + 1) + a(-) J+1/2 (J + |)], (2) 

'-V-' 

Coriolis (K = 1/2) 

where a is the Coriolis decoupling parameter. Electromagnetic transitions among 
the bandmembers likewise follow a well-defined pattern. For the electric quadrupole 
operator, in particular, the reduced matrix element [13] between initial and final band 
members with angular momenta ./, and Jf, respectively, follows the relation 

(^Arlicy***) = (2J* +1) 1/2 (JiK20\J f K) (3) 

v -v-' V -v-' 

Rotational Intrinsic (oc eQ o) 

The value depends on the particular band members involved only through a Clebsch- 
Gordan coefficient, while the details of the intrinsic state determine the intrinsic 
quadrupole moment eQo = (16tt/ 5) Quadrupole moments Q(J), 

proportional to (?/ , ja'||Q 2II'0JA')> and reduced transition probabilities /if M2; J, —y 
Jf), proportional to | (Jf || Q 2 1| J%) | 2 , are uniquely related by (3), to within a common 
normalization determined by Qq. 

4. ROTATION IN THE Be ISOTOPES 

Energies following a rotational pattern are most easily recognized if plotted 
against an angular momentum axis which is scaled as J(J + 1), so that energies 
in an ideal rotational band lie on a straight line (or staggered about a straight line, 
for K = 1/2). Such a plot of energies against angular momentum is shown for the 
calculated states of 9 Be in Fig. 2, both for natural (negative) painty states (left) and 
unnatural (positive) parity states (right). (The parity of the lowest allowed oscillator 
configuration, or traditional shell model valence space, may be termed the natural 
parity, and that obtained by promoting one nucleon by one shell the unnatural par¬ 
ity.) 

These calculations, from Ref. [10], arc obtained taking a specific basis trun¬ 
cation lV max = 10 (for natural parity, or N max = 11 for unnatural parity) and ba¬ 
sis parameter Ml = 22.5 MeV (near the variational minimum). They thus may be 
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Fig. 2 - Energy eigenvalues obtained for states in the natural (left) and unnatural (right) parity spaces 
for 9 Be. Energies are plotted with respect to an angular momentum axis which is scaled to be linear 
in J(J + 1). Solid symbols indicate candidate yrast band members, and shaded symbols indicate 
candidate excited band members. Solid lines indicate the corresponding fits for rotational energies (2) 
(based on the first three band members for K = 1/2 or all band members for K = 3/2). Vertical 
dashed lines indicate the maximal angular momentum accessible within the lowest harmonic oscillator 
configuration (or valence space). 

thought of as taking a snapshot of the spectrum along the path to convergence. 
These calculations are based on the realistic JISP16 nucleon-nucleon interaction [14]. 
The Coulomb interaction has been omitted, to ensure exact conservation of isospin, 
thereby simplifying the spectrum (we consider states of minimal isospin T = T z 
only). However, the effect of the Coulomb interaction, if included, is mainly to in¬ 
duce an overall shift in the binding energies. 

Rotational bands are most readily identifiable near the yrast line, where the 
density of states remains comparatively low. Identification of band members is based 
not only on recognizing rotational energy patterns, but also on observing collective 
enhancement of electric quadrupole transition strengths among band members and 
verifying rotational patterns of electromagnetic moments and transitions. In the nat¬ 
ural parity space [Fig. 2(a)], a K = 3/2 yrast band (solid symbols) and K = 1/2 
excited band (shaded symbols) are identified, while, in the unnatural parity space 
[Fig. 2(b)], a K = 1/2 yrast band is identified (solid symbols). For each of the candi¬ 
date bands, the line indicates a rotational energy fit to the calculated band members. 

Intriguing features of the bands found in the Be isotopes relate to their termi¬ 
nation at finite angular momentum. Some bands terminate at (or below) the maximal 
angular momentum permitted within the valence space, as in the case of both natural 
parity bands of 9 Be [Fig. 2(a)]. [The maximal angular momentum permitted within 
the valence space (or lowest harmonic oscillator configuration) of each parity is in¬ 
dicated by the vertical dashed lines in Fig. 2.] However, some bands also extend 
beyond this angular momentum, as in the case of the unnatural parity yrast band of 
9 Be [Fig. 2(b)]. While this band deviates from the rotational energy formula above 
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Fig. 3 - Dependence of calculated observables on the basis truncation N max , for rotational band mem¬ 
bers in 9 Be: yrast and excited band energies, as (a) absolute energies E and (c) excitation energies E x , 
and yrast band electric quadrupole observables, both (b) quadrupole moments and (d) in-band transi¬ 
tion matrix elements. Values are shown for successive AWx truncations (6 < AI ma x < 10) to provide 
an indication of convergence (calculations of increasing N max are indicated by plusses, crosses, and 
circles, respectively). Electric quadrupole moments and matrix elements are normalized to the intrinsic 
quadrupole moment Q o, as determined from Q( 3/2). The curves indicate the values expected from the 
rotational formula. In panel (d), the upper and lower curves are for A J = 2 and A J = 1 transitions, 
respectively (for the definition and significance of the signs of these matrix elements, see Ref. [10]). 


the maximal valence angular momentum, other calculated bands in the Be isotopes 
continue through this angular momentum with no apparent deviation (see Ref. [10] 
for examples). 

Let us return to the challenge of convergence, now in relation to rotational pat¬ 
terns. The -/V max dependence of the energy eigenvalues, at fixed Ml = 22.5 MeV, is 
shown for the members of both natural parity bands in Fig. 3(a). For each step in 
AM ax- the calculated energies shift lower by several MeV. Flowever, it may also be 
seen that the energies of different members of the same band move downward in ap¬ 
proximate synchrony as N max increases. Thus, the energies within the band remain 
comparatively unchanged, as is seen more directly when we consider excitation ener¬ 
gies in Fig. 3(c). The excitation energy of the excited band relative to the yrast band, 
though not converged, varies much less rapidly with N max than do the eigenvalues 
themselves. 

The strengths of the various electric quadrupole transitions originating from the 
candidate band members in the natural parity bands of 9 Be are shown in Fig. 4. We 
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Fig. 4 - Calculated electric quadrupole transition matrix elements among states in the 9 Be natural 
parity space: (a) from yrast band members and (b) from excited band members. In each panel, shaded 
circles indicate the initial levels being considered for transitions. All angular momentum decreasing 
transitions from the selected levels are shown. Line thicknesses are proportional to the magnitude of 
the reduced matrix element for the transition (also conveyed through a gray scale). Matrix elements are 
shown for the proton quadrupole tensor, but results with the neutron quadrupole tensor are similar. 


observe both the dominance of transitions within each band and the nonnegligible 
cross-talk between bands, as well as to some states outside these bands. 

Although the calculated electric quadrupole observables are far from converged 
[Fig. 1(b)], rotational structure per se is reflected in the ratios of matrix elements 
within a band, not their absolute magnitudes, which depend also upon the intrin¬ 
sic structure via Qo [see (3)]. In Fig. 3 (at right), we consider the calculated elec¬ 
tric quadrupole matrix elements within the 9 Be yrast band — quadrupole moments 
[Fig. 3(b)] and quadrupole transition matrix elements (both A J = 1 and A J = 2) 
[Fig. 3(d)] — and compare with the rotational predictions (lines). The overall nor¬ 
malization Qo is eliminated, by normalizing to the quadrupole moment of the band 
head [from which a value of Qo is determined by (3)]. Thus, rather than consid¬ 
ering the unconverged matrix elements individually, as in Fig. 1(b), we effectively 
only consider ratios, of one unconverged matrix element to another. These normal¬ 
ized matrix elements are seen to be comparatively independent of iV max , at the iV max 
values considered. 

The resemblance between the calculated quadrupole moments and transition 
matrix elements and the expected rotational values (Fig. l(b,d)), while clearly not 
perfect, is strongly suggestive of rotation. One should bear in mind that quadrupole 
moments of arbitrarily chosen states in the spectrum fluctuate not only in magnitude 
but also in sign, and that calculated transition matrix elements for arbitrarily chosen 
pairs of states fluctuate by many orders of magnitude [Fig. 4]. The rotational rela¬ 
tion (3) is equally valid whether we take the quadrupole operator to be the proton 
quadrupole tensor (i.e., the physical electric quadrupole operator) Q ;) or the neutron 
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quadrupole tensor Q n . The matrix elements of these two operators, shown together in 
Fig. 3, provide valuable complementary information for investigating rotation. More 
extensive examples, including bands which extend to higher angular - momentum, and 
considering magnetic dipole observables, may be found in Ref. [10]. 

Returning to the initial questions (from Sec. 1), now that we have explored 
how recognizable the signatures of rotation are, let us quantitatively examine how 
robust the rotational predictions are and how realistic they are in comparison to ex¬ 
periment. The energy parameters for bands across the Be isotopic chain are summa¬ 
rized in Fig. 5: the band excitation energy E x (defined relative to the yrast band as 
E x = E{) — Eo t yrastX th e band rotational parameter or slope A, and the band Coriolis 
decoupling parameter or staggering a (for K = 1/2). Results are shown for a se¬ 
quence of -/V max truncations. Parameters for experimentally observed bands [15] are 
also shown (horizontal lines). (The calculations and experimental data are detailed 
in Ref. [10].) 

Our results for the rotational parameter A and Coriolis decoupling coefficient 
a [Fig. 5(a-c)] are sufficiently stable with respect to lV max to permit at least a rough 
comparison with experiment. The value of the rotational parameter [Fig. 5(a,b)] 
varies by a factor of ~ 2 across the different calculated bands. The calculations 
reproduce both the general range of experimental values for A and the general trend 
of which bands have higher or lower slopes (excepting 1 J Be, where the discrepancies 
may arise due to mixing of the yrast and excited bands [16]). The Coriolis staggering 
for the calculated K = 1/2 bands [Fig. 5(c)] varies in both amplitude and sign, and 
the experimental trend in both these properties is reproduced across the bands. 

However, the excitation energies E x [Fig. 5(d,e)] converge at rates which vary 
considerably among the bands. In general, the excitation energies are decreasing 
with 7V max , bringing them toward the experimental values. Yet, they are varying 
too strongly with iV max for it to be immediately obvious how close the converged 
predictions will lie to the experimental values. 

More detailed comparisons with the experimentally identified rotational bands 
could be made possible by the application of basis extrapolation methods [17-19], 
which are being developed to deduce the converged values for energies (and other 
observables). Such methods are still in their formative stages. Nonetheless, it is 
intriguing to apply a straightforward scheme based on a presumed exponential con¬ 
vergence of energy eigenvalues with IV max [17]. The calculated energies are taken to 
approach the converged value E^, for N max —> oo, as 

E(N maXl hEl) — Ey- Ta/,oexp( bfifiN max ). (4) 

Calculations of the energy at three successive iV max values, for fixed Ml, are sufficient 
to determine all three parameters (af,n, l>i,n, and E 0c ) in (4) and thus provide an 
extrapolation to the converged energy. 



Emergence of rotational collectivity in ab initio no-core configuration interaction cal... 


9 



0.8 

0.6 

> 

u 

°- 4 § 

0.2 


0.0 


30 


20 > 
o 




10 


0 


Fig. 5 - Band energy parameters for odd-mass Be isotopes (left) and even-mass Be isotopes (right): 
the rotational constant A (top), Coriolis decoupling parameter a (middle), and band excitation energy 
E x (bottom). Bands are distinguished as natural parity yrast (circles, solid), natural parity excited 
(circles, shaded), and unnatural parity yrast (diamonds). Values are shown for 6 < AWx < 10 for 
natural parity or 7 < iVmax <11 for unnatural parity, with larger symbols for higher A’ ma x values. 
Parameter values are also shown based on exponentially extrapolated level energies (paired triangles) 
and experimental bands (horizontal lines). 


Using extrapolated energies for the band members, we again determine the 
band energy parameters (paired triangles in Fig. 5), yielding results which match ex¬ 
periment to a greater ( e.g ., 9 Be excitation energies) or lesser ( e.g ., 10 Be excitation 
energy) degree. However, such extrapolations are subject to considerable uncertain¬ 
ties [17]. It is therefore not yet clear to what extent the remaining discrepancies 
reflect deficiencies in the ab initio description of the nucleus, with the chosen inter¬ 
action (JISP16), or limitations of the exponential extrapolation. 

In conclusion, we have seen that, despite the challenges of obtaining converged 
energies or electromagnetic observables, rotational patterns can robustly emerge in 
ab initio NCCI calculations, and these calculations can yield predictions of rotational 
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properties for both qualitative and quantitative comparison with experiment. What 
remains to be addressed, from among our initial questions (in Sec. 1), is the physical 
origin and intrinsic structure of the rotation: prospective contributors include multi¬ 
shell SU(3) symmetry (shown to be dominant in NCCI calculations [6]), symplectic 
symmetry [3,20], and a-a molecular rotation [21]. 
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